Excited states in poly-diacetylene chains: 
A Density-matrix-renormalization-group study 
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We study theoretically poly-diacetylene chains diluted in their monomer matrix. We employ the 
density-matrix renormalization group method (DMRG) on finite chains to calculate the ground 
state and low-lying excitations of the corresponding Peierls-Hubbard-Ohno Hamiltonian which is 
characterized by the electron transfer amplitude to between nearest neighbors, by the electron- 
phonon coupling constant a, by the Hubbard interaction U, and by the long-range interaction V. We 
treat the lattice relaxation in the adiabatic limit, i.e., we calculate the polaronic lattice distortions for 
each excited state. Using chains with up to 102 lattice sites, we can safely perform the extrapolation 
to the thermodynamic limit for the ground-state energy and conformation, the single-particle gap, 
and the energies of the singlet exciton, the triplet ground state, and the optical excitation of the 
triplet ground state. The corresponding gaps are known with high precision from experiment. We 
determine a coherent parameter set (tl = 2.4 eV, a* = 3.4eV/A, U* = 6eV,V* = 3eV) from a fit 
of the experimental gap energies to the theoretical values which we obtain for 81 parameter points 
in the four dimensional search space (to, a, U, V). We identify dark in-gap states in the singlet and 
triplet sectors as seen in experiment. Using a fairly stiff spring constant, the length of our unit cell 
is about one percent larger than its experimental value. 

PACS numbers: 71.20.Rv, 71.10.Fd, 78.30.Jw, 78.20.Bh 



I. INTRODUCTION 

Poly-diacetylene (PDA) chains dispersed with low 
concentration in their monomer single-crystal matrix 
are prototypical quasi one-dimensional materials. i"— The 
structural disorder in the chains and their surrounding 
matrix is tiny, the electronic excitation energies of the 
diacetylene monomers are much higher than those of the 
polymer, and the chains' electronic excitations in the en- 
ergy range of visible light can be measured with a very 
high accuracyi^ 

Exciton-polaritons have been generated that have been 
shown to be coherent over tens of micrometers, i.e., 
several ten thousand monomer unitsi^ This observation 
was confirmed by weight measurements after dissolving 
the chains and their monomer matrix.'^ Consequently, 
the opto-electronic properties of the PDAs result from 
the electrons' mutual interaction and their interaction 
with the lattice potential, while the influence of disor- 
der is negligible. This makes these materials the perfect 
testing-ground for theoretical model studies which de- 
scribe interacting electrons on perfectly ordered chains. 

The typical single-particle gap in PDAs is Eg^p ^ 
2.4 eV, see Sect.|TTl Density-functional theory band struc- 
ture calculations in the local-density approximation for 
generic PDA geometries estimate the bare band-gap to be 
-E^barogap ~ 1.2 cV, Or lessi^ Rcsults from various meth- 
ods are compiled in table 1 of Ref. fi\; recent calcula- 
tions using the Perdew-Burke-Ernzerhof global hybrid 
density functional and the 6-311G(2d,2p) basis set of 
atom-centered Gaussian functions (geometries from the 



TPSS density functional and the 6-31G(o?) basis set)^ 
give -Bbaregap ~ 1.6 cV.^ The comparison shows that elec- 
tronic exchange and correlations account for a substan- 
tial fraction of the single-particle gap. In contrast to 
inorganic semiconductors, the exciton binding energy in 
PDAs amounts to about 20% of the single-particle gap. 
Such large binding energies suggest that the electron- 
electron interaction must be treated accurately for the 
calculation of the optical properties of the PDAs. 

In order to describe the optical excitations in PDAs, 
two approaches have been taken. The first approach 
starts from an ab-initio density-functional theory calcu- 
lation of the bare band structure in local-density ap- 
proximation (LDA), which is then supplemented by an 
approximate treatment of the residual electron-electron 
interaction, e.g., the GW approximation for the single- 
particle bands and the Bethe-Salpeter equation (BSE) for 
the excitons {LJ:iA+GW+BSE),i2Jl Actual calculations 
for the PDAs often omit the GW step (Wannier the- 
ory).!- Within this approach, a number of experimental 
data can be reproduced, e.g., the l^B^ exciton binding 
energy and its polarizability. 

This approach is less successful for the triplet sector 
Typically, the energy of the triplet ground state is too 
high. Recall that, if the electron-electron interaction is 
absent, the energy of the triplet ground state is iden- 
tical to the single-particle gap. Starting from a weak- 
coupling description of the electron-electron interaction 
on the chains, it is difhcult to obtain the experimentally 
observed energy renormalization by a factor of almost 
three. Moreover, polaronic effects are not considered in 
the LDA-I-BSE approach. 
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The second approach to a theoretical description of 
the primary excitations in polymers starts from a many- 
particle model Hamiltonian that describes only the tt- 
electrons and their mutual interaction. Typically, empir- 
ical parameters are used for the tight-binding band struc- 
ture and for the Pariser-Parr-Pople (PPP) potential;^ 
With the help of the density-matrix renormalization- 
group (DMRG) methodjiS- the ground state and elemen- 
tary excitations for such models can be calculated for 
large chains with very high accuracy. In this way, the 
electron-electron interaction is treated without resorting 
to any approximations. 

In a recent study^^ we used the Hubbard-Ohno poten- 
tial and the tight-binding parameters of Ref. to cal- 
culate the binding energy, polarizability and wave func- 
tion of the singlet exciton, in good agreement with ex- 
periment. However, in our previous study we could not 
reproduce satisfactorily the energy of the triplet ground 
state. Moreover, we did not take polaronic effects into 
account. Here, we shall overcome these shortcomings. 

In this work, we perform an extensive DMRG study 
of the Peierls-Hubbard-Ohno Hamiltonian for the tt- 
electrons on a chain. We start from the tight-binding 
Peierls description of Race et al but replace their Ohno 
potential^ by the Hubbard-Ohno potential.— The essen- 
tial difference between the two parameterizations of the 
Pariser-Parr-Pople interaction^'' lies in the treatment of 
the Coulomb interaction for 7r-electrons on the carbon 
atom. The local Hubbard repulsion potential is substan- 
tially larger than the corresponding Ohno interaction.— 

Since we are mostly interested in the polaronic ef- 
fects, we ignore the energetic effects introduced by the 
ligands R,R'. Some preliminary studies show that the 
energy difference between the singlet and triplet ground 
states is not influenced by the introduction of a local 
potential for the carbon atoms which are linked to the 
side groups. We presume that the dominant influence of 
the side groups comes from the presence (or absence) of 
strain in the chains. As can be seen from the data for 
3BCMU and 4BCMU, different ligands and the resuU- 
ing strain results in small energetic differences, of about 
0.1 eV for the nBCMU family. 

The outline of our work is as follows. In Sect. [Hi 
we summarize the experimental observations on the sin- 
glet and triplet in-gap states in PDAs. In Sect. IIIH we 
define the Peierls-Hubbard-Ohno Hamiltonian and the 
model parameters which provide the basis of our numer- 
ical DMRG study. In Sect. |TVl we briefly discuss our 
numerical approach. In Sect. |Vl we motivate the param- 
eter regime that we choose for our study. In Sect. |Vll we 
present our results. In Sect. IVlIl we summarize and con- 
clude. Technical details are deferred to the appendices. 



II. EXPERIMENTAL OBSERVATIONS 

We start with an overview of the experimental obser- 
vations relevant for our study. 



A. Ground-state conformation 

First, we collect relevant experimental data on the 
ground-state properties. 

1. Lewis structure 

The diacetylene monomer building unit is comprised of 
four carbon atoms. The four outer electrons of each car- 
bon atom are hybridized. Three of them form localized 
bonds. There are cr-bonds between neighboring carbon 
atoms on the chain. Two carbon atoms on the chain share 
a local TT-bond made by two Py electrons. The other two 
carbon atoms share cr-bonds to covalent ligands R and 
R' , which are several Angstr0m long and differ for vari- 
ous members of the PDA family. In this work, we focus 
on poly-(Butoxy-Carbonyl-nMethylene-Urethane) (poly- 
nBCMU) chains with n = 3, 4 where the side groups are 
given hy R = R' = (CH2)„ - OCONH - CH2 - COO - 

(CH2)3CH3. 

The fourth carbon electron is delocalized over the car- 
bon backbone in a molecular 7r-orbital. Due to the Peierls 
effect, in the ground state the tt electrons dimerize the 
chain into an alternating sequence of short and long 
bonds. After dimerization, the four carbon atoms in the 
unit cell are linked by a triple bond, a single bond, a dou- 
ble bond, and a single bond. The corresponding Lewis 
structure of the ground state is shown in Fig. [TJ 

r R 1 

L R' Jn 

FIG. 1. Lewis structure of a poly-diacetylene unit cell. 



2. Lattice parameters 

For a high-quality single crystal of Poly-[l,2-bis-(p- 
tolylsulfonyloxymethylen)-l-buten-3-inylen] (PTS), the 
atomic distances at room temperature have been mea- 
sured^l as n = 1.191(4) A, Td = 1.356(4) A, and = 
1.428(4) A [uncertainties in the last digit in brackets] for 
the triple (t), double (d), and single (s) bonds, respec- 
tively. Typical atomic distances for other PDA poly- 
mer single crystals are2^ Vt = 1.20 A, — 1.36 A, and 
rs = 1.43 A. The same set of data applies for 3BCMU- 
PDA at low temperatures f22, The chain of atoms is not 
straight; the single and double bonds alternately form an- 
gles of (fii = 120° and ip2 = 240° degrees, with a tempera- 
ture variation of a few degreesi^ In the comparison with 
our calculations we shall assume that the bond lengths 
and angles given above for the PDA single crystals are 
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representative for chains in their diacetylene monomer 
matrix. 

We shall only deal with planar ('blue') PDA chains.^ 
The individual polymer chains of 4BCMU are strained 
in their monomer single crystals but are essentially un- 
strained in SBCMUp^ Strain should be the primary 
source for differences in the spectra of these two PDAs. 



B. Excited states 

Next, we summarize experimental results on the low- 
lying electronic excitations. PDAs are center-symmetric 
insulators. Their ground state G is a spin singlet with 
symmetry Ag under inversion. 



1. Single-particle gap and excited singlet states 

The charge gap for single-particle excitations, as de- 
termined from Franz-Keldysh oscillations in electro- 
absorption experiments, 2" is £;gap(3BCMU) = 2.482 eV 
and £:gap(4BCMU) = 2.378 eV in 3BCMU and 4BCMU 
chains, respectively. 

The excitation energy of the primary singlet exciton S 
(symmetry ^B^) defines the optical gap, A^p^ — Eg — 
Eg, which amounts to A^pt(3BCMU) = 1.896 eV and 
A^pt(4BCMU) = 1.810 eV in 3BCMU and 4BCMU, re- 
spectively. Therefore, the singlet exciton binding energy, 
defined by A|^ = Eg^p - Es, becomes A^^(3BCMU) = 
0.586 eV in 3BCMU, and A^^(4BCMU) = 0.568 eV in 
4BCMU, about 24% of the band-gap. The energy levels 
are sketched in Fig. [21 

The singlet exciton is the energetically lowest state in 
the spin-singlet sector which can be generated by a single- 
photon absorption process. In addition, there are further, 
optically dark states in the gap. The existence of opti- 
cally dark states Xi and X2 below the optical gap can 
be inferred from non-radiative decay processes which are 
monitored via pump-probe spectroscopy, see Ref. 3] for 
a review. The exciton rapidly populates the states Xi^2 
so that they should have the same spin quantum number. 

In principle, the energy of singlet states in the gap 
with ^^g-symmetry can be determined via two-photon 
absorption. Two-photon absorption for a single-crystal of 
the poly-diacetylene paratoluene-sulfonate reveals three 
gap states with energies Ex^ = 1.05Aqpj, Ex2 = O.OA^p^, 
and Exi = OSA^^^. These states should exist in all 
PDAs. Note that in the experimental literature;^ the 
numbering of the states Xi and X2 is reversed. 



S. Phonon energies 

Raman scattering reveals vibrational energies which 
are assigned to the oscillations of the double (D) and 
triple (T) bonds. For 3BCMU chains in their monomer 
matrix they are ?iWD = 0.181 eV and Hlut = 0.261 eV.^'^ 
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FIG. 2. Energy levels of in-gap states in the spin-singlet 
and spin-triplet sectors. Single-tip arrows: optical absorption 
spectroscopy; double-tip arrows: two-photon absorption spec- 
troscopy. Double arrows: binding energies (gaps). G: singlet 
ground state (l^Ag); S: singlet exciton (l^Bu); Xi, X2, X3: 
singlet dark states (m^Ag); T: triplet ground state (l^Bu); 
T*: optical excitation of the triplet ground state (I'^Ag); Y: 



dark triplet state (m By)- 



In accordance with the Raman data, the optical spectra 
of 3BCMU chains show strong exciton replicas at the en- 
ergies Ec — Es+huju — 2.079 eV and at Ed = Es+huj^ — 
2.160 eV, respectively, when the exciton is accompanied 
by single optical phonons corresponding to the vibrations 
with frequencies wd and wt- 

Electro-absorption measurements^!^ show that there 
are more significant single-phonon replicas of the singlet 
exciton at the energies Ea.t = Es + Hujs.b' with ?ia;s = 
0.090 eV and Tiojb- = 0.155 eV. Due to the fact that 
Eb ~ Ec, the phonon replica at Eh appears in the low- 
energy flank of the replica at energy Ec- The electro- 
absorption measurements also permit the identification 
of multi-phonon replicas, e.g., at replica energies ?i(ti;D* + 
wd), 2huJY), h{uJu +Wt), 2huJT^ 



3. Triplet ground state and excited triplet states 

The triplet sector is more difficult to access experimen- 
tally because a transition between the spin-zero ground 
state G and the lowest spin-one state T is optically for- 
bidden. Optical pump-probe spectroscopy^ reveals that 
a small fraction of singlet excitons decays into a long- 
lived state which can be optically excited by the probe 
pulse. Its long life-time indicates that this in-gap state 
is the spin-triplet ground state T. The probe pulse gen- 
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erates transitions from T to T* in the triplet sector. The 
optical gap in the triplet sector is defined as A^p^ = 
Et- ~ Et, and amounts to A^pt(3BCMU) = 1.360 eV 
and A*pt(4BCMU) = 1.345 eV, respectively. 

Optical pumping above a threshold Ef « 2.0 eV very 
efficiently generates states which show a strong optical 
absorption with energy A^p^. This can be readily un- 
derstood if the spin-singlet excitations above E{ fission 
into triplet pairs. In turn, these triplet pairs can re- 
combine into singlets and decay optically. If we ignore 
lattice effects (bi-polaron formation), we obtain a rea- 
sonable estimate for the energy of the triplet state T, 
A,t = Et - Eg Pi Ef/2. For 3BCMU and 4BCMU 
this estimate gives Ast(3BCMU) > 1.0 ± 0.05 eV and 
A,t(4BCMU) > 0.95 ± 0.05 eV, respectively. 

The binding energy of the triplet ground state is de- 
fined by A*^ = Eg^p-Er- It amounts to A^^(3BCMU) « 
1.5 eV and A*^(4BCMU) « 1.4 eV in 3BCMU and 
4BCMU, respectively, more than 60% of the single- 
particle gap. The energy of the optically excited triplet 
ground state T* is found to be Et* (3BCMU) - Eg = 
2.36 eV and £;t.(4BCMU) ~ Eg ^ 2.30 eV above the 
ground state, about 0.1 eV below the threshold i?gap for 
single-particle excitations. 

The optically dark singlet states Xi and X2 decay non- 
radiatively into a state Y. Its weak population and long 
life-time indicate that it is reached via intersystem cross- 
ing so that it ought to be a spin-triplet state which lies 
energetically above the triplet ground state T. It should 
be of symmetry so that it cannot be reached via an 
optical excitation of the triplet ground state. 

Fig. [2] shows the experimentally observed level spec- 
trum for PDA chains in the singlet and triplet sectors. 
We summarize the corresponding values for the in-gap 
states in table HI and compare them to our theoretical 
results for our best parameter set, (ig = 2.4 cV, a* = 
3.4eV/A, U* = 6eY,V* = 3eV), see Sect. ED 



III. MODEL DESCRIPTION OF 
POLY-DIACETYLENE CHAINS 

In this work we restrict ourselves to the description 
of the TT electrons because they dominate the optical re- 
sponse of the poly-diacetylene chains immersed in their 
monomer matrix for energies hw < 3eV. In order to 
make contact with previous work,^^ we treat the other 
electrons as inert, i.e., they are supposed to form the un- 
relaxed geometry of the carbon backbone. The distance 
between two carbon atoms is R2 for a single cr-bond, and 
Ri < i?2 for the a-py double bond (extrinsic distortion). 



A. Electronic Hamiltonian 

The motion of 7r-electrons between neighboring carbon 
atoms and their mutual Coulomb interactions defines the 



Enerj 






3BCMU 


4BCMU 


Theory 


Ex I 






1.5 


1.4 


1.74 [1.94] 








1.7 


1.6 


1.85 [1.94]'^ 


Es = 


Aopt 




1.896 


1.810 


2.00 [2.05] 


Ex3 






2.0 


1.9 










2.482 


2.378 


2.45 [2.47] 




- E 




0.586 


0.568 


n 45 42 


Et 


Ast 




1.0 ± 0.05 


0.95 ± 0.05 


1.00 [1.06] 


/ = 


= Ast - 




2.36 ± 0.05 


2.30 ± 0.05 


2.25 


A* 




1.360 


1.345 


1.25 [1.28] 


A* - 


- -E^gap 


- Et 


1.5 ± 0.05 


1.4 ± 0.05 


1.45 [1.40] 



^ Several degenerate states are found in DMRG. 



TABLE I. First and second column: Excitation energies in 
3BCMU and 4BCMU at low temperatures. All energies are 
measured in eV relative to the energy of the ground state, 
-Eg = 0. Bold number: directly measured; italic number: es- 
timate. Third column: our results, see Sect. IVll the numbers 
in square brackets give the excitation energy for the rigid- 
lattice transition from G (-Egap, -Es, E-Xi 2, ^t) and from T 
(AU). 

electronic problem, 

Ho=f + V, (1) 

where T and V specify the electrons' kinetic energy and 
their mutual interaction, respectively. 

1. Ktnetic energy 

The motion of the 7r-clectrons over the unrelaxed back- 
bone is described by the operator for the kinetic energy, 

T^~^tl {claCl+l,a + ci^i^^h^a) , (2) 
l-cr 

where cf^, ci,a are creation and annihilation operators, 
respectively, for a 7r-electron with spin a =t,i on site I 
with two-dimensional coordinate n = (xijUi)"^. The ma- 
trix elements ti are the electron transfer amplitudes be- 
tween neighboring sites. Transfer amplitudes between 
next-nearest neighbors should be included to fit better 
the band structure of all carbon electrons.^ The ampli- 
tude for an electron transfer between two carbon sites at 
distance tq = 1.4 A is given by which we use as an 
adjustable parameter. 

We consider the half- filled band exclusively, i.e., in the 
ground state and for the excitations in Fig. [2] the number 
of TT electrons Ne equals the number of lattice sites N. 

2. Coulomb interaction 

The diacetylene monomer single-crystals are insula- 
tors, and also the PDA chains display a finite charge gap. 
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Therefore, the long-range Coulomb interaction is not dy- 
namically screened at the energy scale of a few electron 
volts. 

Therefore, we start from the Pariser-Parr-Pople (PPP) 
interaction^^ 



V = U 



N 

E 

1=1 



N 

E 



[{fli - 1) {flrn - 1)] 



(3) 



Here, hi — nj^-j- + n/,^ counts the number of electrons on 
site I, and hi^a — cl„ci,ij is the local density operator at 
site / for spin a. The strength of the (local) Hubbard 
interaction is parameterized by U , and Vi^m are the PPP 
parameters for the effective Coulomb repulsion between 
electrons at different positions n and f„i- 

For the description of electrons and holes in quantum 
wires and other quasi one-dimensional structures in vac- 
uum, various effective potentials have been used in the 
literature i^^'^^ For example, in our previous studyi^ we 
used the erf-potential {x =\ri — fm\) 

Vi% = y-f(x) = ^^exp[{x/Rf] [1 - erf(x/i?)] , 

(4) 

where erf(a;) is the error function, R is the adjustable 
confinement parameter, and = 2.3 is the static di- 
electric constant for the diacetylene monomer matrix. In 
general, the PPP interaction for the 7r-electrons on the 
chain in the surrounding matrix has the form 

Ed 



^ppp ^ ^ppp 



(x) = V"'{x) 



e{x) 



(5) 



where e{x) is the static dielectric function at distance 
X = In — rVnl with e{x oo) — e^. Unfortunately, the 
short-distance behavior of e{x) is unknown. 

In this work we follow Refs. [H, [13, and H^] and ap- 
proximate the Pariser-Parr-Pople interaction using the 
Hubbard-Ohno potential, i.e., for x = \fi — rm\ ^ we 
set 

V 



V, 



PPP 

Lm 



^Qhno ^ yOh„o(^) ^ 



V 



14.397 eV 



i + i3{x/ky 



(6) 



The Ohno potential and its adjustable parameter V de- 
scribe the effective strength of the Coulomb interaction 
at short distances; for large electron-electron distances, 
yOhno(^ ^ oo) ^ / (cdx) bccausc = 14.397 eV A. 
Note that U and V are independent adjustable parame- 
ters in our theory. Later, we shall assume that the screen- 
ing of the on-site interaction is substantially less effective 
than for the long-range interaction. 

The eigenstates of the electronic problem follow from 
the solution of the corresponding Schrodinger equation. 



(7) 



We denote expectation values of operators A in the nor- 
malized state I*) as (A) = 



3. Particle-hole symmetry 

The Hamiltonian is invariant under the particle- 
hole transformation q^^ (— l)'c/^_o.- half band- 
filling, the ground state |G) is also invariant under this 
transformation. 

In our numerical investigation, we calculate important 
excited states using the proper quantum numbers for the 
spin symmetry {S = 0,1; p = 2S + 1), inversion symme- 
try {X = Ag,Bu), and particle-hole symmetry {v = ±1). 
Therefore, we label states in the form m ^X"" where 
m > 1 counts the states with the same symmetry in as- 
cending energetic order. The quantum numbers for the 
most prominent in-gap states are summarized in table HIl 





State 


Symmetry 


G 


T 


S 


T* 


Spin (5) 





1 





1 


Inversion [X) 


1 


-1 


-1 


1 


Particle-hole (v) 


1 


1 


-1 


-1 


Classification 




1 







TABLE II. Quantum numbers of important in-gap states. For 
the definition of the states, see Fig. [21 



B. Electron-lattice interaction 

The electron-phonon coupling leads to the dimeriza- 
tion of the ground-state structure (Peierls effect). More- 
over, excitations carry a polaron cloud, and the polaronic 
shifts in the single-particle excitation energies will in gen- 
eral be different from those of bound pairs. For example, 
the singlet exciton fissions into a bound pair of triplet 
polarons. This bipolaron has an energy which is lower 
than twice the energy of a triplet polaron. Therefore, 
the estimates in table |T] provide lower bounds for Agt . 

Due to the (small) Peierls distortion, the energy in- 
creases so that the total Hamiltonian reads 



H = H^ 



N-l 

E 

1=1 



(8) 



At distance tq = 1.4 A, the electron transfer matrix el- 
ement is given by our adjustable parameter to. The 
strength of the electron-lattice coupling is parameterized 
by the coupling constan t^^'^^ 



A, 



2a' 



(9) 



where a is the strength of the Peierls coupling and Ki are 
the elastic constants for the carbon backbone, namely. 
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K for the cr-bonds and G > K for the a-py bonds, Ki = 
K + (5; mod 4, lie — K). The parameters a, K, and G 
must be adjusted. In effect, we address the region where 
A; < 0.1 so that the adiabatic approximation is valid, i.e., 
we may treat the lattice distortions classically. 

The intrinsic Peierls distortion implies a modulation of 
the bond lengths. 



5ri = ri - Ri = - 



2a 



(10) 



Here, i?2 ~ rg for a single cr-bond of the carbon backbone 
chain and Ri = tq — 5^ / (2a) for its cr-py double bond. 
The size of the extrinsic dimerization 5'^ is calculated in 
App. [X] As a result of the intrinsic Peierls dimerization, 
the ground-state unit cell of the distorted chain consists 
of four carbon atoms, linked by a single bond, a double 
bond, a single bond, and a triple bond, see Fig. [1] 

The energies 5i modulate the electron transfer ampli- 
tudes. 



ti = ti + 5i/2 



(11) 



with 



ti=to + (572 

(undistorted (J-py double bond), (12) 

U = to 

(undistorted a single bond). (13) 

Note that, for small distortions, the transfer matrix el- 
ements for the single, double, and triple bonds in the 
ground state obey 



(14) 



where Tg, r^, and are the lengths of the single, double, 
and triple bond in the unit cell, see eq. ([T0|) . 



C. Model parameters 

Our model employs the following parameters: (i) The 
electron transfer matrix element to for a single C-C- 
bond at distance tq = 1.4 A; (ii) the strength of the 
local Hubbard interaction U; (iii) the strength of the 
short-range Coulomb interaction V; (iv) the Peierls cou- 
pling a; (v) the spring constants K and G. The model 
parameters are adjusted to reproduce the single-particle 
gap -Egap , the singlet exciton energy Eg , the energy of the 
triplet state E^, and the optical gap in the triplet sec- 
tor Aqpj, see table m Moreover, we estimate the values 
for the spring constants from the energy of the phonon 
replicas. 

Of course, we cannot scan a five-dimensional param- 
eter space completely. Therefore, we have to restrict 
ourselves to values which seem plausible, see Sect. El 
The model parameters investigated are summarized in 
table Hm 



Fixed 
parameter 


Value 


Control 
parameter 


Range of 
values 


Ed 

K 

G 
ro 


2.3 

44 eV/A^ 
68 eV/A^ 
1.4 A 


to 
V 
U 
a 


2.0 eV . . . 2.4 eV 
2 eV . . . 3 eV 
5 . . . 6 eV 
3.4 eV/A ... 3.6 eV/A 



TABLE III. Parameters of the Peierls-Hubbard-Ohno model. 



IV. METHOD 

First, we outline our procedure to find the optimal lat- 
tice structure. Next, we define the single-particle gap and 
describe how we address excited in-gap states. Lastly, we 
remark on our DMRG procedure and the extrapolation 
of our finite-size data to the thermodynamic limit. 

In this section and in the remainder of the paper, all 
energies {to, U, V; E'gap, A^p^, A^t, A'pJ are given in eV, 
all lengths are given in A, and a is given in units of eV/A. 



A. Optimization of the lattice structure 

1. Procedure 

The values for the electron transfer amplitude modu- 
lations follow from minimization of the energy functional 
E^i,{5i) = {^\H\^) for the normalized state l^") which can 
be the ground state or any excited state of H. The actual 
values for the dimerization are obtained from the mini- 
mization of the energy functional Exj,(6i) with respect to 
6i subject to the constraint 



(15) 



1=1 



This reflects the fact that the total chain length should be 
fixed; see below. The condition is taken into account 
with the help of the Lagrange multiplier F. 

According to the Hellmann-Feynman theorem)^ the 
negative derivatives of the energy functional with respect 
to Sri define the force fields //, 



1 ^ 



/31//A2 



2ed 



dx 



^[l-|-/3(|rl-r;|/A)2]3/ 



(17) 



y^'' ds, 



((n, - 1) (h, - 1)) 



where 6p {p = 1, . . . , N — 1) are the Peierls modulations 
of the electron transfer amplitudes 

yi,3 =yi- Vj- 
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For fixed bond angles, xi and yi are defined as xi = 
xi-i+ri cos(60°) and yi — yi-i — ri sin(60°) for tlie bonds 
at an angle of 120°, and xi = xi^i + ri and yi — yi-i 
otherwise, with ri = Ri — Si/{2a) from ([TU]). The force 
fields fi are zero at the optimal values S°^^ for a chosen 
state l^"). Note that \'i>) is an eigenstate of the electronic 
problem which is parameterized in terms of Sp. There- 
fore, the minimization of the force fields has to be done 
self-consistentlyJ^ 

i. In step k of the iteration {k — 1,2, . . .), the target 
eigenstate l^'fe), e.g., Gk, Sk, or Tk, is calculated 
for Sk-i using the infinite-lattice DMRG algorithm. 
In all our cases, the initial choice 6i-i = for Z = 
1, . . . , — 1 leads to converged solutions. 

ii. For given k and fixed quantum-mechanical expec- 
tation values in the distortion energies are 
determined iteratively. 

To this end, the condition fi — in (ITBl) is used 
to determine the distortion energies for the next 
iteration, S]^+^ = 27r<oA,(-r^ [^^.p]), (n > 0, 
^k-i — ^k;i)- Here, the Lagrange parameter follows 
from (HSi as = Y^i ^iPi[Slp\/ Ei ^i- 

The distortion energies typically converge after 
some five to fifteen iterations. The converged solu- 
tion defines Sk+i-i = lim„_i.oo (5^.; for the next iter- 
ation in k. 

iii. The steps (i) and (ii) are repeated until a converged 
set of distortion energies and DMRG energies for 
the states are obtained. Si = lirak^oo Sk-i, |^) = 
limfc_^oo 

Strictly speaking, the condition of a fixed chain length 
ic corresponds to a^i jv + 2/i at — ^1- We have verified 
numerically that the condition (|15p preserves the chain 
length up to 0.1% for N < 100. 

2. Polaronic energies 

It is important to note that we optimize the lattice 
structure and the corresponding distortion energies for 
each state separately. Our excited states contain all po- 
laronic energy contributions, i.e., we give their relaxed 
energies. This polaronic relaxation was not taken into 
account in our previous studyJ^ There, we studied rigid- 
lattice transitions with fixed electron transfer amplitudes 
ts, td, and tt, which correspond to the Lewis structure of 
Fig. [TJ In general, they are higher in energy than the 
corresponding relaxed excitations. 

It is not a priori clear whether the relaxed or the rigid- 
lattice energies should be compared to experiment. For 
the optical singlet excitation, no Stokes shift is observed 
between absorption and fluorescence spectre^ so that, in 
the Franck-Condon picture, the exciton creation process 
corresponds to a vertical transition. 



In order to estimate the polaronic contribution to the 
energy, we calculate the energy of excited states in the 
rigid-lattice approximation for our optimal parameter 
set, see Sect. El For {t^ = 2.4 eV, a* = 3.4eV/A,f/* = 
6 eV, V* — 3 eV), the relaxed energy of the single-particle 
gap is £'g°'p'^'"* = 2.45 eV whereas the energy of a single- 
particle excitation with flxed electron transfer matrix ele- 
ments tf, t'^, and t^ leads to E'^g'^ = 2.47 eV, see tableHl 
Thus, the energy relaxation due to the polaron formation 
amounts to about ^poiaron = 0.02 eV. In the band pic- 
ture, the singlet exciton is a bound state of particle-hole 
excitations. Correspondingly, the polaronic shift in Eg 
should be about twice as large as <5poiaron, as indeed ob- 
served, i;j,<=ia'^<=d ^ 2.00 eV and E^'^''^ = 2.05 eV, so that 

^rigid _ ^relaxed ^ q.OS cV « 2<5polaron. 

The same amount of polaronic relaxation energy is 
observed for the singlet-triplet gap, iJ^'®"^ — i^^ciaxcd _ 
0.06 eV. When we start from the relaxed triplet ground 
state, we find for the optically excited state T* that 
E'^f'^ = 1.28 eV whereas E'^^}'''""^ = 1.25 eV. It is seen 
that the polaronic relaxation energy amounts to about 
(^polaron also in the triplet sector. 

Our observation of a fairly small polaronic relaxation 
energy ties in with the fact that our electron-lattice cou- 
pling is small. A; < 0.1, and the adiabatic approximation 
is valid. In our comparison with experiment below, we 
show the energies for transitions between lattice-relaxed 
configurations. 

B. Single- particle gap and in-gap excitations 

The band-gap or single-particle gap -Egap is defined by 
the difference in chemical potentials for a system with A'e 
and Nc ~ 1 particles, 

E^gap = KNc) - KNc - 1) , 
H{N,) ^ Eg{N, + 1) - Eg{N,) , (18) 

where Eq{Nc) is the energy of the A'e-particle ground 
state G. In the presence of particle-hole symmetry at half 
band-filling, we have 

£;gap = 2fi{N,) (19) 

for the minimal energy of a single-particle excitation. 

In poly-diacetylenes, the singlet exciton and its vi- 
bronic replicas carry most of the oscillator strength of 
the optical excitations. The quadratic Stark effect in the 
electro-absorption proves that they are bound states of 
electron-hole excitationsi^ The exciton energy thus de- 
fines the optical gap, 

A^pt = Es {N, ^ N) - Eg{N, ^ N) , (20) 

where Eg {Nc — N) is the energy of the first excited state 
of the half-filled system with symmetry i?„ . The binding 
energy of the exciton is then obtained as 

= i?gap - Es . (21) 



8 



Note that we calculate for finite-size systems so that all 
quantities must be extrapolated into the thermodynamic 
limit, N oo. 

For a full account of all in-gap states, we target up 
to five states simultaneously in the spin-singlet and spin- 
triplet sectors, respectively. Note that the lattice relax- 
ation must be done for each state separately. These cal- 
culations represent the most time consuming part of our 
investigations. 



C. Numerical procedure 

In this work we present results from numerical density- 
matrix renormalization-group (DMRG)i^ calculations on 
finite chains with open boundary condition (OBC) us- 
ing an adopted version used previously J-^ The discarded 
weight was kept below 77 = 10~^ for all calculations by 
employing the dynamical block-state selection (DBSS) 
procedurcj^i^^ We have set the minimum number of 
block states to Mmin = 400 and used three sweeps. The 
maximum value of the number of block states varied 
around Afniax = 600. As benchmarks we compared our 
DMRG energies for some selected parameter values with 
those from a DMRG code used earlier by Race, Bar- 
ford, and Bursill.^^ The latter one, however, only uses 
the infinite-lattice procedure so that our variational en- 
ergies are always slightly lower. 

The ground state as well as all excited states have been 
targeted and relaxed individually using fc = 5 to 15 re- 
laxation iteration steps to reach the pre-set convergence 
criterion on di. Note that each relaxation iteration step 
requires a full DMRG run with three sweeps. The whole 
relaxation procedure has been performed for all target 
states and for all chain lengths independently from = 6 
up to iV = 66, in steps of AiV — 4. For our optimized 
parameter set we have performed calculations for up to 
N = 102 sites. 

The spin multiplicity of the converged target states 
was calculated from the expectation value of the operator 
for the total spin S*^, employing the expectation values 
of the corresponding correlation functions.'^'* In order to 
determine the optically dark in-gap singlet states we have 
shifted the energies of the triplet states out of the gap by 
adding the term '^■^ '^^ '^j" ^'^ Hamiltonian. As an 
alternative procedure, we have also identified the exciton 
states by calculating the dipole strength as introduced 
in Ref. [1^ in which the reduced density matrix of the 
target state was constructed from the reduced density 
matrices of the ten lowest eigenstates. 

For the 81 parameter points of the four-dimensional 
search space (to, a, U, V) and for the energies of the states 
(G, S, ^Jgap, T, T*,Xi,Yi) shown in Fig. H we have taken 
on average A: = 10 relaxation steps for the 16 to 24 dif- 
ferent chain lengths for single and multiple target states 
which results in about 130.000 fuU DMRG runs. We 
estimate that the calculations consumed overall about 
45 CPU years which were provided by 100 parallelized 




0.04 0.06 
1/N 



FIG. 3. (Color online) Finite-size scaling of the charge gap 
Bgap, calculated for to = 2.4 eV, a = 3.4 eV /A, U = 6eV, 
V = 2.0, 2.5, 3.0 eV, and 14 < Af < 66 sites. The lines are 
quadratic fits. 




0.04 0.06 
1/N 



FIG. 4. (Color online) Finite-size scaling of the singlet exciton 
energy Es (optical gap A^pt), calculated for to — 2.4 eV, a = 
3.4eV/A, (7 = 6eV, y = 2.0, 2.5, 3.0 eV, and 10 < iV < 66 
sites. The lines are quadratic fits. 



CPUs so that the calculations took some five months in 
real time. 

The PDAs are charge and spin insulators, i.e., the gaps 
for single-particle, optical, and magnetic excitations are 
finite. The materials are characterized by finite correla- 
tion lengths. Therefore, end effects decay exponentially, 
and local operators that are calculated in the middle of 
the chain display a regular behavior as a function of in- 
verse system size. Thus, various quantities that we cal- 
culate for finite chain lengths A^ can be extrapolated re- 
liably to the thermodynamic limit, A^ — >■ cxi, by using a 
second-order polynomial fit. 

As an example, in Fig. [3] we show the charge gap cal- 
culated for to = 2.4, a = 3.4eV/A, U = 6eV, and 
V = 2.0, 2.5, 3.0 eV, as a function of 1/A^. As seen from 
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the figure, the second-order polynomial fits permit a reli- 
able extrapolation but the quadratic curvature becomes 
dominant for TV > 50 only. Therefore, a study of long 
chains is mandatory. The appearance of the inflection 
point at sizable chain lengths is more pronounced for the 
gap states than for the single-particle gap. An example, 
the optical gap, is shown in Fig. 01 




o V=2 
> V=2.5 
° V=3 



■ 0.02 0.04 0.06 0.08 0.1 

1/N 

FIG. 5. (Color online) Finite-size scaling of the triplet ground- 
state energy Et (singlet-triplet gap Ast), calculated for to — 
2.4 eV, a = 3.4eV/A, !7 = 6eV, = 2.0, 2.5, 3.0 eV, and 
10 < A'' < 66 sites. The lines are quadratic fits. 




o V=2 
> V=2.5 



■ 0.02 0.04 0.06 0.08 0.1 

1/N 

FIG. 6. (Color online) Finite-size scaling of the optical gap in 
the triplet sector A*pt — Et* —Et, calculated for to ~ 2.4 eV, 
a = 3.4 eV/A, U = 6eY,V = 2.0, 2.5, 3.0 eV, and 14 < TV < 
66 sites. The lines are quadratic fits. 

In Figs. [5] and |6] we show finite-size results for the 
triplet ground-state energy Et (singlet-triplet gap Agt) 
and for the optical gap in the triplet sector A*p^ = 
Et* - Et, respectively, for to = 2.4 eV, a = 3.4eV/A, 
U = 6eV, V = 2.0, 2.5, 3.0 eV. The ground-state energy 
Et in the triplet sector rapidly converges as a function 



of inverse system size 1/N because the state is deep in 
the gap. The energy Et* of its optical excitation T* is 
close to the threshold Sgap for single-particle excitations 
so that long chains must be studied for a reliable ex- 
trapolation to the thermodynamic limit. For the scan of 
our parameter regime as specified in the next section, we 
limit ourselves to chains of length A'^ < 66. The accuracy 
of the extrapolation is better than SE = 0.05 eV. 



3 



° V=3 






= V=3 



0.05 001 0.05 0.1 
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FIG. 7. (Color online) Finite-size scaling of (a) the charge gap 
Egap, (b) the optical gap A^pt, (c) the singlet triplet gap Ast, 
and (d) the optical gap in the triplet sector Aopt, calculated 
for the optimal parameter set {tg — 2.4 eV, a* = 3.4eV/A, 
U* = 6eV, V = 3eV), and N < 102 sites. The lines are 
quadratic fits. The extrapolated values are summarized in 
Tab.lll 

Before we detail our optimization of the parameter set 
in Sect. IVI Al we show in Fig. [7] the finite-size scaling 
of all four gaps for our optimal choice of parameters 
{t*a = 2.4 eV, a* = 3.4eV/A,[/* = 6eY,V* = 3eV) 
for which we investigate chains with up to A^ = 102 
sites. As can be seen from Fig. [71 the finite-size ex- 
trapolation can be done very accurately with an uncer- 
tainty of SE ^ 0.01 eV. The extrapolated values are 
Eg^p = 2.45 eV, A^p^ = 2.00 eV, A,t = 1.00 eV, and 
Aj^pj = 1.25 eV, see also tableHl 

V. PARAMETER REGIME 

Before we present our results in Sect. IVH we give ar- 
guments for the parameter regime used in our study, as 
summarized in table IIIIl 

A. Spring constants 

At first sight, it seems to be easy to obtain the spring 
constants K and G because they can be inferred from Ra- 
man scattering data for short molecules. Unfortunately, 
the influence of the delocalized tt electrons appears to 
be crucial. For example, early Raman experiments in 
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ethane^ H3C-CH3, give = SleV/A^ whereas for 
the single-bond in diacetylene^aS HC=C-C=CH, /^^^ = 
45eV/A^ is found. In polymers, the situation is equally 
ambiguous because the theoretical analysis of the same 
Raman data for poly-acetylene leads to the same set of 
concurring values^ K^^-^ = 31eV/A^ and if^^'^ = 
46eV/A^i^ In the present work, we investigate the con- 
sequences of a strong spring constant for the cr-bond, 
K > 40eV/A^. We plan to present a detailed study of 
the vibrational properties in the near future. 

In order to estimate K (and G), we calculate the opti- 
cal phonon spectrum from a simple classical model. The 
carbon atoms of mass M = 12u in the Lewis structure 
of Fig. [1] are linked by spring constants of strength K 
and G — gK in the sequence (G, K, K, K) in the unit cell, 
corresponding to the a-py-bond and the three <7-bonds. 
The optical phonons for two-dimensional vibrations of 
the chain are derived in App. |B] There are four positive 
solutions of the characteristic equation for the phonon 
frequencies, < = yj2K/M < ujc < ujd, where uJb 
is the resonance frequency of two carbon atoms linked 
by the spring constant K.2S~ The frequencies uja^c,d are 
obtained from the zeros of the third-order polynomial 

p{y) = _ 2{g + 2)y^ + (7/2 + 6g)y ~ 3g (22) 



as uJa.cA = ^/vZIdKjM, p{ya,c,d) = 0. 

In comparison with experiment, see Sect. Ill we as- 
sign ujd{g) = and ujcig) = wd- From (clIt/wd)^ = 
(0.261/0.181)2 = 2.079 we find go = 1.547, yd{go) = 
4.463, ydgo) = 2.147, and ya{go) = 0.4845. From 
ydgo) = MojI,/K we obtain K = 44.1 eV/A^ which 
agrees with results obtained for poly-acetylene chains^^ 
^PA,2 ^ 46eV/A2. In addition, we find G = gK = 
68.3 eV/A^ for the a-py-hond. For comparison, the 
spring constant in ethene (ethylene) was derived as 

The other two optical phonons have the energies huJb — 
n^/WjM = 0.175 eV, and fiuja = n^K/M^ya{go) = 
0.086 eV, respectively. These values are in good agree- 
ment with experiment, TiUb ~ Tiwd* — 0.155 eV, and 
hiva ~ hujs = 0.090 eV, see Sect. IIIB21 Note that we 
use experimental data for comparison which include the 
influence of the 7r-electrons whereas for our model cal- 
culations we employ bare values for the backbone. The 
influence of the itinerant tt electrons must be calculated 
self-consistently so that the values for K and G need fur- 
ther refinement. This task is left for a future study. 



B. Electron-phonon coupling 

Next, we discuss the bare bandstructure for non- 
interacting electrons and estimate the size of the electron- 
phonon coupling constant a. 



1. Bare bandstructure 

The bare bandstructure for the ground state with 
filled valence bands with energies E^^2{k) = — e2(fc) and 
Ev.i{k) = — ei(fc) and empty conduction bands with en- 
ergies Ecs{k) = ei(fc) and Ec,2{k) = f-2{k) is derived in 
App. O see eq. (jUgj) . The bare gap is given by A'^''™ = 
2ei(fc = 0). The Coulomb interaction enhances all gapsi^ 
so that the singlet-triplet gap Agt — Et — Eq w 1 eV 
will be larger than A^^'^°. Therefore, the size of the bare 
band-gap constrains the possible values for the electron- 
phonon coupling. 

In the absence of Coulomb interactions, the values for 
is, td, and tt must be determined from the minimization 
of the total ground-state energy per unit cell, 

etot('5d,(5t) = ekin('5d, (5f) +epot{Sd,St) , 

where g = G/K and A — 20^ /{ntoK). The kinetic en- 
ergy of the electrons is given by 

— ie,ik)+e2ik)), (24) 

where the factor two accounts for the spin degeneracy. 

The electron transfer matrix elements tg.d.t and the 
distortion corrections 6s^d,t are related by tg — to — [St + 
<^d)/4, td = to + Sd/2, and' it = <o + St/2 + 6^2. Here, we 
used the fact that 2ds + Sd + 6t = because the length 
of the unit cell is not changed by the intrinsic distortion. 
The strength of the extrinsic dimerization (5° follows from 
the solution of eq. (KE\i of App. |X] for U = V = 0, 

6%U = V = 0) _ 8Ti(<5°) _ 8T,{6'^) _ 

TTtoX 2eo{S-) 4Ti(<5-) ^ ' 

so that S%U = F = 0) Aa^/K. 

The numerical minimization of etot(^d, (^t), eq. (|23l) . 
leads to the somewhat surprising result that the bare 
gap A^^'^°{to,a) very weakly depends on t^. For a = 
3.5eV/A, K = 44eV/A2, and G = 68eV/A2, we 
find that A'^^''<'(2,3.5) = 0.878 eV and A'^'^^<=(2.4, 3.5) = 
0.835 eV; it even decreases slightly with increasing to. 
In contrast, the bare band-gap strongly increases as a 
function of a. For to = 2.4 eV, K = MeV/k"^, and 
G = 68eV/A2, we find A'^'^''<'(2.4, 3.4) = 0.776 eV and 
^barc(2.4^3.6) = 0.897 eV. The bare band-gap becomes 
larger than Ast = 1 eV for to = 2.4 eV and a = 3.8eV/A. 
Therefore, we must use smaller values for a as derived 
and used previously] ^^i^^ 

2. Intrinsic and extrinsic Peierls distortion 

In poly-acetylene (PA), the mobile 7r-electrons dimer- 
ize the chain. This intrinsic Peierls effect results in a mea- 
sured bond length alternation of Ar = 0.04 A, i.e., long 
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FIG. 8. (Color online) Regions in the parameter space to-V 
which lead to Ri = 1.25 A for U = 6eV. The color coding 
gives the appropriate value for a. 



and short bonds of length rf^ — 1.44 A and r^^ = 1.36 A 
alternate along the chaini^ Almost the same amount of 
alternation in and is seen in PDAs, see Sect. Ill A 21 
Previous DMRG studies^* lead to Ar = 0.03 A. 

When we assume that the intrinsic Peierls effect af- 
fects the triple bond in the same way as the double 
bond, we come to the conclusion that the length Ri 
of the cr-py-bond before the intrinsic dimerization is 
i?i « 1.25 A assuming rt = 1.21 A. Therefore, the 
extrinsic dimerization due to the p^-bond accounts for 
Ta — Ri = '5'^/(2a) = 0.15 A. For non-interacting elec- 
trons, we have 5°{U — V = Q) — Aa^/K, independent 
of to- Thus, we arrive at the estimate a ~ 0.15A(iir/2) 
which gives a « 3.3eV/A for K = 44eV/A2. 

Even in the presence of Coulomb interactions, the cal- 
culation of 5° is simple because it requires the solution of 
a two-site problem only, see App. [Xj Therefore, a com- 
plete parameter scan is readily accomplished. As we 
shall argue below, the on-site Coulomb repulsion U is 
quite substantial. Therefore, in Fig. |S] we show param- 
eter regions in {tQ,a,V) for U = 6eV that correspond 
to S°/{2a) = 0.15 A. As compared to non-interacting 
electrons, the electron-phonon coupling a has to be in- 
creased by some 10% to generate the same extrinsic 
Peierls dimerization for the interacting two-site system. 
This indicates that the Coulomb interaction makes the 
bonds noticeably stiffer. In order to account for this ef- 
fect, our spring constant K from Sect. IV Al should be 
reduced by at least ten percent; a more thorough scan 
for the if-parameter will be done in a future study. 

Fig.[5]shows that the dependence on to and V is rather 
weak: it is mostly the Peierls coupling a that determines 
the size of the bond length shift. From the data in Fig. [8] 
we conclude that a « (3.5 ± 0.1)eV/A is a reasonable 
starting point, whereby we compensate our somewhat 
too large spring constant K. 



C. Coulomb parameters 

Previous studiesiSiiiii^ii^ succeeded to reproduce the 
charge gap and the energy of the singlet exciton. How- 
ever, the complexity of the in-gap states could not be 
recovered. In particular, the splitting Agt of the sin- 
glet and triplet ground states cannot be reproduced as it 
comes out substantially too high. Moreover, in previous 
approaches no dark singlet states have been found that 
lie energetically below the singlet exciton. 

The energetic positions of the states T and Xi_2 are a 
clear signal of substantial electronic correlations induced 
by the Hubbard interaction U. Our initial calculations 
with small ratios U/V put T and Xi 2 too high in en- 
ergy as compared to experiment 4S As we shall show in 
the next section, we find a reasonably good description of 
the level scheme in Fig.[2]only for substantial U and com- 
parably smaU y , 5 eV < [/ < 6 eV and 2 eV < V" < 3 eV, 
as indicated in table IIIII so that U/V = k k, 2 holds. 
Note that we included the dielectric constant explic- 
itly in the Ohno potential ^ because we treat chains 
immersed in their monomer matrix. 

Substantial values for the Coulomb interaction were 
advanced by Chandross and Mazumdar>22 as a result of 
their model study of poly-phenylene-vinylene (PPV) thin 
films. In order to describe the linear and non-linear op- 
tical properties of PPV, they proposed Uqm — 8eV with 
C^Cm/Vcm = k = 2i^i^ Note, however, that these au- 
thors worked with fixed lattice parameters, i.e., without 
lattice relaxations for the excitations, and employed an 
approximation (single configuration interaction, SCI) to 
calculate optical excitations. 

The bare bandwidth in our calculations \sW k, Ai^ > 
9eV which still is larger than the on-site interaction. 
Therefore, the system is still far from the spin-Peierls 
limit. At the same time, however, the correlations are 
strong enough to impede weak-coupling approaches;^ 



VI. RESULTS 

First, we scan our parameter space and determine our 
best parameter set (to,a, C/, F). Next, we analyze the 
energy levels of optically dark in-gap states and comment 
on the lattice parameters. 



A. Optimization of the parameter set 

For each choice of the parameter set (to,a, C/, F), we 
calculate the ground-state energies at half band-filling 
and one additional particle, and the energy of the three 
excited state S, T, and T* at half band-filling from which 
we determine the four gaps i?gap , A^p^ , Agt , and A* pj for 
systems with size 10 < TV < 66. The lattice geometry of 
all states is relaxed, see Sect. IIVI We performed calcula- 
tions for (^0, Q!, [/, V) for a broad range of parameters us- 
ing small system sizes up to iV = 30 sites. For the optimal 
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range (to = 2.0, 2.2, 2.3, 2.4 eV, a = 3.4, 3.5, 3.6 eV/A, 
U = 5.0, 5.5, 6.0 eV, V = 2.0, 2.5, 3.0 eV) we investigated 
systems with up to = 66 sites, i.e., we address alto- 
gether 108 different parameter sets. 

As an example, in Fig. |9] we show the extrapolated 
energies for fixed (to — 2.4 eV, U = 6.0 eV), as a func- 
tion oi V — 2.0, 2.5, 3.0 eV for the three parameters 
a = 3.4, 3.5, 3.6 eV/A. As expected for our small pa- 
rameter window, we observe a fairly linear dependence 
of the gaps on the parameters V and a. It is seen that 
not all the gaps can be reproduced perfectly with a single 
parameter set. In general, for fixed {to,U), the optical 
gap in the triplet sector, A*pt = Et' —Et requires larger 
values for {a, V) than the other gaps. Therefore, we have 
to compromise to find a good parameter set. 




V 



FIG. 9. (Color online) Extrapolated energies as a function 
of y for a = 3.4,3.5, 3.6 eV/A for fixed (to = 2AeV,U = 
6.0 eV). The horizontal lines give the experimental results for 
3BCMU (green dashed) and 4BCMU (black dotted). 

To this end, we define V* as the value V for a 
given (to, a, U) which reproduces the experimental gaps i 
{i = 1,2,3,4 for Eg^p, Al^„ A,t, A^pJ for 3BCMU and 
4BCMU from table HI Then, we calculate the joint stan- 
dard deviation CT(to, a, U) from 

[a(to,a,[/)]^ ^}-Y}{V*~V*)/eYY . (26) 

The optimal parameter set minimizes the spread cr. The 
results for a(to,ck,f/) are shown in table HV] for 3BCMU 
and for 4BCMU in table |Vl respectively. 

The two tables IIVI and |V] indicate that the best set 
for both 3BCMU and 4BCMU is (to = 2.4 eV, a = 
3.4eV/A, t7 = 6eV). Then, a look at Fig. [9] shows that 
V = ?> eV is the best value for which we have data avail- 
able. Therefore, we shall use the set (tg = 2.4 eV, a* — 
3.4eV/A,C/* = 6eV,l/* = 3eV) as our optimal param- 
eter set. The trend shows that a might even be a bit 
smaller, a < 3.4eV/A, and to a bit larger, to ^ 2.4 eV, 
for the optimal case. We plan to perform a more system- 



U = 5.0 eV 


to/eV 
2.0 2.2 2.4 


a = 3.4eV/A 
a = 3.5eV/A 
a = 3.6eV/A 

U = 5.5 eV 


1.96 2.79 3.22 
1.68 3.12 3.67 
2.66 3.63 4.11 

to/eV 
2.0 2.2 2.4 


a = 3.4eV/A 
Q = 3.5eV/A 
a = 3.6eV/A 

U = 6.0 eV 


1.91 1.63 2.28 
1.96 2.14 2.85 
2.53 2.63 3.37 

to/eV 
2.0 2.2 2.4 


a = 3.4eV/A 
a = 3.5eV/A 
a = 3.6eV/A 


14.07 1.95 1.55 
3.56 2.06 2.11 
3.59 2.91 3.67 



TABLE IV. Spread cr as a function of (to, Q, U) for 3BCMU 
where the experimentally observed gaps are given by -Egap = 
2.482 eV, A=pt = 1.896 eV, A^t = 1.0 eV, and A|;pt = 
1.360 eV. The three best sets (to, a, U) are printed in bold. 



u 


= 5.0 eV 


to/eV 
2.0 2.2 


2.4 


a = 


3.4eV/A 


2.42 3.19 


3.66 


a — 


3.5eV/A 


2.30 3.65 


4.16 


a = 


3.6eV/A 


3.36 4.17 4.58 


U 


= 5.5 eV 


to/eV 
2.0 2.2 


2.4 


a — 


3.4eV/A 


2.85 2.31 


2.91 


a — 


3.5eV/A 


2.67 2.87 3.45 


a = 


3.6eV/A 


3.37 3.34 


3.96 


U 


= 6.0 eV 


to/eV 
2.0 2.2 


2.4 


a = 


3.4eV/A 


17.89 2.51 


2.15 


a = 


3.5eV/A 


3.96 2.56 


2.68 


a — 


3.6eV/A 


4.07 3.64 


5.12 



TABLE V. Spread as a function of (to, a, t/) for 4BCMU 
where the experimentally observed gaps are given by -Egap = 
2.378 eV, A=pt = 1.810 eV, A^t = 0.95 eV, and A^pt = 
1.345 eV. The three best sets {to,a,U) are printed in bold. 

atic estimate for the optimal parameter set in the near 
future. 

For our optimal parameter set (tp = 2.4 eV, a* = 
3.4 eV/A, U* =6 eV, V* = 3 eV), the corresponding the- 
oretical values for the excitation energies are -Egap = 
2.45 eV, A^pt = 2.00 eV, A^t = 1.00 eV, and A^^^ = 
1.25 eV. They are also given in table HI 



B. Optically dark in-gap triplet states 

In the triplet sector, we find a series of optically dark 
states Y; {I = 1, 2, 3, 4) at energies just above the triplet 
ground state. As indicated in Fig. ^ the optically dark 
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FIG. 10. (Color online) Scaled excitation energies of the 
first four triplet states, Di = {N/A + l)(i?Y, — Et), as a 
function of inverse system size 1/A'^ for 10 < TV < 78 and 
{t*o = 2.4 eV, a* = 3.4 eV, (7* = 6eY,V* = 3eV). The 
crossings of the ordinate are equidistant, Di = 0.637 meV, 
D2 = 1.235 meV, D3 = 1.851 meV, and D4 = 2.548 meV, 
indicating a linear dispersion relation for the excitations. 



triplet states Y; open decay channels for the states Xi 
and X2. 

The structure of the optically dark in-gap triplet states 
is also interesting from a theoretical point of view. In a 
band picture, the triplet ground state with spin compo- 
nent = 1 consists of a hole in the | valence band and 
an electron in the f conduction band, both at momentum 
fc = 0. Close in energy are the corresponding excitations 
at finite but small momentum k. The dispersion of these 
excitations is quadratic as a function of k. 

In the interacting case, these excitations may form a 
spin-flip density wave with momentum q whose disper- 
sion relation at low-energy is given by 

Csf ((?) = CsiQ , (27) 

where Cgf is the (sound) velocity. For finite chains, we 
have quantized quasi-momenta, qi — Trl/[{N/4 + l)d] {I — 
1,2,...), where d is the length of the unit cell. Therefore, 
the levels Y; should obey 

EY,~ET^eM + j^ + ... , (28) 
(j + l)(£^Y,-i?T) = C3f5 + ^ + ... . (29) 

We confirm this hypothesis in Fig. [10] where we show 
{N/4: + 1)(-By, - Et) as a function of 1/N. As seen 
from the figure, the energy differences scale to a finite 
value linearly in 1/A^. Moreover, the extrapolated val- 
ues are equidistant from which we can read off Agf = 
Csfir/d = 0.625 meV. Using d = 4.9 A we thus esti- 
mate the velocity for the spin-flip density excitations 
as Cst = 148 m/s. For localized spin models, we have 
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FIG. 11. (Color online) (a) Energies of the lowest singlet 
states as a function of inverse system size 1/A^ for 10 < A'' < 
78 for {t*o = 2.4 eV, a* = 3.4eV,f/* = 6eV,V* = 3eV). 
After the extrapolation, the states X2a and X2ti are essen- 
tially degenerate, Exi = 1.744 eV, i?x2a = 1.853 eV, Ex2b = 
1.863 eV, and -Es = 1.995 eV. The horizontal dashed line gives 
the value for the charge gap in 3BCMU. (b) Overlap intensity 
of the dipole operator for spin-singlet in-gap states after the 
extrinsic and the intrinsic relaxation for A'' = 46 sites. 

cj « J(d/4) so that, approximating cj ~ Csf, the effec- 
tive magnetic interaction in our system is of the order of 
one meV, J « (4/7r)Asf = 0.8 meV. 



C. Optically dark in-gap singlet states 

In Fig. [TTk we show the (relaxed) excitation ener- 
gies for the four energetically lowest singlet states as a 
function of 1/iV for {t* = 2.4 eV, a* = 3.4eV/A,[/* = 
6eV,V* = 3eV). For our optimal parameter set, there 
are at least three optically dark singlet states below the 
singlet exciton. We extrapolate i?Xi = 1.744 eV which 
is about 0.25 eV higher than estimated from experiment, 
see tablelH This indicates that the local correlations could 
still be larger than U = 6eV. Note, however, that the 
absolute positions of Xi and X2 have not been deter- 
mined experimentally for nBCMU but it is known from 
pump-probe spectroscopy'^ that there are (at least) two 
optically dark singlet states below the singlet exciton. 

The next two dark in-gap singlets are almost degener- 
ate in energy, i^Xaa = 1.853 eV and i^Xab = l-863eV. 
They lie below the singlet exciton, as seen in experi- 
ment, but about 0.15 eV higher than estimated in tableU] 
Fig. [Tib shows the dipole overlapi^ for the four singlet 
states Xi, X2Q, X26, and S at chain length N = 46. Only 
the exciton S has a flnite dipole overlap. Note the energy 
shift of all in-gap states due to the intrinsic Peierls effect. 
The Peierls shift amounts to several tenths of an eV. 

As in the triplet sector, we employ the band picture 
at flxed particle number to view the elementary excita- 
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tions of the ground state as a hole in the valence band 
and an electron in the conduction band. In a Wannier 
picture, the electron-electron interaction forms bound 
states from these pairs, such as Xi, the lowest- lying Ag 
singlet, and S, the singlet exciton with i?„ symmetry. 
Apart from these bound states, there should be a contin- 
uum of scattering states. The degeneracy of the states 
and X2{, indicates that they are near the thresh- 
old to the X-continuum. Indeed, two-photon absorption 
above Ef « 2.0 eV excites states that can fission into two 
triplets/^ Unfortunately, it takes a significant amount of 
CPU time in our DMRG approach to target more than 
four in-gap singlet states simultaneously so that a more 
detailed investigation of the in-gap spectrum remains an 
open problem. 



D. Lattice parameters 

Finally, we show results for the lattice constants in 
Fig. [T^l For our optimal set of parameters, the ex- 
trapolated values are rg = 1.425 A, rd = 1.373 A, and 
rt = 1.239 A, for the single, double, and triple bond, 
respectively. When compared to the experimental val- 
ues for PDA single crystals given in Sect. Ill Al the val- 
ues for the single and double bonds are rather good but 
the triple bond is too large. This can also be seen from 
the size of the unit cell. Given the Lewis structure of 
Fig. [TJ the unit cell in chain direction has the length d 



with (f = rl + (2rs -|- n)^ - 2rd(2rs -I- rt) cos{tpi). Using 
ifi = 120°, this results in d = 4.92 A which is slightly 
larger than the experimental values for 3BCMU chains, 
rfsBCMU = 4.89 A4 



1.35- 




1.25- 



0.08 
1/N 



FIG. 12. (Color online) Atomic distances Vs, Vd and n for 
the single, double, and triple bonds in the chain center as 
function of inverse system size 1/A'^ for 6 < A'^ < 66 for {to = 
2.4eV,a = 3.4eV/A,C/ = 6eV), and V = 2.0, 2.5, 3.0 eV. 
The extrapolated values for the best parameter set = 3 eV) 
are n = 1.425 A, rd = 1.373 A, and n = 1.239 A. 



The comparison shows that the spring constants K 
and G are too large and/or the electron-phonon coupling 
constant a is too small. Consequently, the spring con- 
stants K, G should be included as parameters in the 
optimization procedure. 



VII. CONCLUSIONS 

A. Parameter values 

In our study we use quite sizable values for the 
Coulomb parameters. Substantial values for the Hub- 
bard interaction, Ucm = 8 eV, were used by Chandross, 
Mazumdar et al. in their studies of the optical proper- 
ties of poly-phenylene-vinylene (PPV) thin filmsi^ For 
PDA chains immersed in their monomer matrix, we find 
U — 6 eV and confirm their previous result for the ra- 
tio between U and V, J7cm/Vcm = k = U/V = 2^2dl 
Note that we additionally screen the long-range part of 
the Coulomb interaction by the dielectric constant of the 
monomer matrix, ed = 2.3. A further increase of the 
Hubbard interaction beyond U — 6 eV would be prob- 
lematic, as can be seen from Fig. [T5k . For fixed to and 
a, the spread of the gaps increases with increasing U. In 
particular, the calculated singlet-exciton energy would 
deviate significantly from its experimental value. More- 
over, the number of dark singlet states below the exciton 
would become larger than expected from experiment. 
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13. (Color online) Optimal values V* for the four gaps 



gap, ^opt5 



Ast, and A^pt (a) as a function of U for to — 



2.4eVanda = 3.4eV/A, (b) as a function of a for to = 2.4 eV 
and U = 6eV, and (c) as a function of to for U = 6eV and 
a = 3.4eV/A. 

For fixed to = 2.4 eV and U = 6eV, a change from 
a = 3.4eV/A would increase the spread in the gaps, 
as can be seen from Fig. [T5b . Therefore, we advocate 
an electron-phonon coupling strength which is some 15% 
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smaller than the value proposed by Ehrenfreund et al., 
aE = 4.0eV/Aja 

In the literature, typical values for the electron transfer 
integral are to = 2.5 eV-^^'^ and <o = 2Ae\f^^^. Our 
work indicates that a readjustment of the bare electron 
transfer integral is not necessary for nBCMU. This can 
be seen from Fig. [T3t, where the spread of Vi increases 
for both smaller and larger Iq. 



B. Outlook 

Our study can be extended in several directions. First, 
the strength of the spring constants K and G must be 
determined self-consistently. In this way, a better agree- 
ment between theory and experiment for the lattice pa- 
rameters can be obtained. 

Second, the parameter search in the four-dimensional 
space {tQ,a,U,V) can be optimized by performing a 
smooth interpolation between the 81 data sets which 
we have investigated numerically thus far. Moreover, 
the present analysis does not distinguish between vari- 
ous ligand types (3BCMU, 4BCMU). The differences are 
partly due to strain, and a distinction between PDAs 
must therefore be incorporated in the lattice parameters. 
In addition, different side-groups introduce (small) elec- 
trostatic potential at the carbon atoms which are linked 
to the side groups. Our numerical analysis shows that 
the influence of an electrostatic potential of the order of 
Ci < 0.3 eV at the side-group sites changes the spectra by 
6E ^ 0.05 eV so that this contribution can be ignored to 
first approximation. 

Third, for an optimal parameter set it is interesting to 
study the properties of the ground and excited states in 
more detail. For example, the lattice structure of the ex- 
citon (polaron-exciton) and its polarizability can be cal- 
culatedi^i^ and be compared with experiment. More- 
over, further in-gap states should be address, e.g., the 
spin-2 ground state l^Ag, or a second singlet exciton 

These tasks are left for a future study. 
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Appendix A: Extrinsic dimerization 

Following Ref. [l^ , we diagonalize the two-site Peierls- 
Hubbard-Ohno model in the spin singlet sector to derive 
the ground-state energy. To this end we consider the two 
states in site representation 



|2> - /I i^tA 



iC2,t 1 I vac) 



i + ^tt^ll) I vac) 



In the subspace 5 = 0, the Hamilton matrix of the 
electronic problem has the entries {Hc)i.j = {i\Hc\j) 
= 1,2). A short calculation gives 



where 



y((5'=) 



-U/2 -2Ti{6'=) 
~2Ti{S'=) C//2-F((5<=) 



(Al) 



i?i(<5^) = ro-^ , Ti(^^)=to + y . (A2) 

i?i((5°) is the length of the cr-Py double bond, Ti((5°) is 
the corresponding electron transfer amplitude, U is the 
strength of the electrons' local Coulomb repulsion, and 
V{5'^) is their interaction on neighboring sites. Note that 
in Ref. [l^, the unshifted energies were used, i.e., the pre- 
vious expressions follow from ours after an energy shift 
by U/2 + V{5°). 

The ground-state energy eo of the two-electron system 
follows from the diagonalization of the matrix ()Aip as 



(A3) 

Now that we know the ground-state energy of the elec- 
tronic Hamiltonian explicitly, we do not have to invoke 
the Hellmann-Feynman theorem in order to determine 
the optimal values for the electron transfer amplitudes (5° . 
Instead, we directly minimize the total ground-state en- 
ergy 



Ea{5-) = eo(<5°) + 



47rioA 



(A4) 



with respect to 5° (A = 2a^/(7r_ftrto)). We set 



1 



2aedA2 [l + ;3(i?i(,5<=)/A)2]3/2 ■ 

(A5) 

Therefore, the optimization of the electron-lattice prob- 
lem for two carbon atoms with a double bond leads to 
the implicit equation 



[2eQ{5'') + U]V'{5^)-m{6'') 



ntoX 2eo((5°) + T/(5<=) 

which is solved iteratively. 



(A6) 
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Appendix B: Optical phonons 



2. Equations of motion 



Assuming the Lewis structure of Fig. [l] the position 
of the atoms in the unit ceh are denoted by the two- 
dimensional vectors Ai, Bi, Ci, and _D/ for / = 1 , . . . , i = 
iV/4. Their equilibrium positions are denoted as A; q, 
Blo, Clq, and Di^. The PDA structure implies B^q — 
Aifl = Difl-Ci,o = rse^, Cifl-Bi^o = ncx, and A;+i.o- 
Di,o — Td{cos{<p)ex — siTi{(j))ey) for the singlet, triplet, and 
doublet bonds, where </> = 180° - (pi = 60°. 



1. Lagrange function 

We shall treat the atomic motions classically. For con- 
venience, we use periodic boundary conditions, L + 1 = 1. 
The kinetic energy of the atoms is given by 



L r 



1=1 L 



(Bl) 



In the spring-constant model, the atoms' potential energy 
is approximated by {g = G/K) 



1=1 



Bi-Ai 



Di-Ci 



Ci-B, 



Ai+i - Di 



(B2) 

To second order in the displacement 5x — x — xq we can 
write 



(1^1 - l^ol)^ ~ (Sx ■ Xa/\xo\ 



(B3) 



We define the (small) displacements ai ^ Ai — Ai^ — 
afe^ + a^Cy, bi = Bi-Bi^ = bfe^ + bfey, q = Ci~Cifl = 
cf etc -I- cfcy, and di = Di - Di^o = dfe^ + dfey. Then, 
the Lagrange function in the harmonic approximation 
becomes L = T — V with 



T = 



1=1 L 



-ft 



1=1 



{bf-aff+g{cf~bff + {df-cff 



+ [cos((/.) - df) - sin(0) (af^, - df)]' 



(B4) 



The variables 5f and cf are cyclic and drop out of the 
problem. 



The Euler-Lagrange equations can be solved using the 
Fourier Ansatz 



xiit) 



k=l 



(B5) 



with k — 27rmfe/L, nik — 0, l,...,i — 1 as the crys- 
tal momentum. Here, ^ = {a^ , ,"f^ , , S^) cor- 
responds to a; = {a^ ,0^ ,b^ ,d^ ,dy). The resulting 
set of six algebraic equations can be cast into a ma- 
trix equation, M{uj)Ek = (0,0,0,0,0,0)^ for the vector 
Ek = {a.%,a\, 131,^^,5^,5^)'^ . With the abbreviations 
y — MuP' jK, c — cos((/)), and s — sin(^), the matrix 
reads 

(B6) 




— cs 



V 



H = 












y ~cs 


-1 







— CS 


~y 










-1 


1- 




-9 











~g l- 




y -1 









-1 


1 1 2 

1 + c- - y 












—cs 



-y 



The equation for the vibrational frequencies w„(fc) as 
a function of the crystal momentum k results from the 
characteristic equation, det(A^(a;)) = 0. 

We are interested in the optical phonon modes, w„ = 
^n{k — 0). The characteristic equation reduces to 



(B7) 



p{y) = y' - 2{g + 2)y^ + (7/2 + 6g)y - 3g 

for all g and (j) = 60°. The finite- frequency solutions are 
denoted as w„ (n = a, b, c, d). We set ojb = \/2K/M, and 
o;^ < o;^ < result from the three real roots of p{y) ~ 
in (HZl) as Wa,c,d = \/~KyZZdW- 

Appendix C: Bare dispersion relation 

The operator for the kinetic energy T in eq. ^ is 
readily diagonalized for periodic boundary conditions, 
iV -I- 1 = 1. For the ground state of non-interacting elec- 
trons, the unit cell consists of four sites, N = AL, and the 
electron transfer amplitudes follow the periodic pattern 
{ts, tt, ts, td). Thus, we may write 

+isC4„+3 ,^C4„+4_o- + i(iC4„^4 ,jC4(„+i) + i^cr^ + h.C. . 

(CI) 

We introduce the four operators bM;k,a for electrons with 
quasi- momentum k = — 7r-|-27rmfe/L, ruk = 0, 1, . . . , L — 1 
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bM: 



1 

1 J2 



The inverse transformation reads 



C4n+M,< 



(C2) 



(C3) 



where we used that the sum over the quasi-momenta k 
generates the orthogonahty relation 



(C4) 



for two lattice indices m, n. In turn, the sum over lattice 
indices n leads to the orthogonality relation 



L-l 



(C5) 



for two quasi-momenta k,p. 

When we apply the transformation (jC3P to the kinetic 
energy (ICip . we obtain 



The remaining task is to band-diagonalize the kinetic en- 
ergy. To this end we diagonalize the 4 x 4-niatrix Aik 
with 



f -t, -tde-"'\ 
-ts -tt 
-tt -t. 



(C7) 



Its eigenvalues in ascending order are Ei{k) = — e2(fc), 
E2{k) = -ei(fc), E:i{k) = ei(fc), Ei{k) = e2(k) with 



ei(fc) < e2(fc). We findil 

[e^,2ik)r^ts+tl/2 + t^j2 



±\J{tl - ^?)V4 + tl [tl + tl + 2tit, cos(fc)] . 

(C8) 



We write 



Mk=U+VkUk 



(C9) 



where V = diag(— e2(A:), — ei(fc), ei(fc), e2(fc)) is a diago- 
nal matrix which contains the eigenvalues of A^fe. We 
define the band-diagonal operators 



Pk,. 

Ik.a 



and find 



/ 6l;/c,cr \ 
fe2;*:,cr 
b3;k,a 



(CIO) 



[K,Jk,a^'^k,a^k,a 



(Cll) 



The a- and /3-bands are the valence bands which are filled 
in the ground state at half band-filling. The 7- and 6- 
bands are the conduction bands which are empty in the 
half-filled ground state. 

The bare gap at half band-filling obeys A'^'*'''= = 2ei(0) 
with 



[61 (0)]2 = i2 + tl/2 + tll2 - VW+K^Witd +tt)/2. 

(C12) 

For small deviations, is,id,is ~ to, this simplifies tc"'^^ 



\tt + ~ 2ts 



(C13) 
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